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ALGEBRAIC  DESCRIPTION  OF  FORMS  IN  SPACE 


by 

Chester  E.  Jenserr^^ 

The  expected  spatial  relation  between  continuous  variables 
is  often  expressed  graphically  early  in  a regression  or  correla- 
tion analysis.  In  turn,  the  graph  may  be  expressed  algebraically. 
To  quantify  or  statistically  evaluate  the  relation,  the  algebraic 
form  may  be  fitted  mathematically  to  a relevant  data  set.  How  to 
approximate  the  graph  algebraically  in  such  a process  is  the  con- 
cern here. 

Given  a dependable  basis  for  the  expected  relation,  curve 
form  bias  introduced  algebraically  and  associated  information  loss 
will  tend  to  diminish  as  the  algebraic  form  approaches  that  of  the 
graph.  When  the  analyst's  confidence  in  the  graphed  form  is  great, 
he  should  then  be  prompted  to  adopt  a high  standard  of  descriptive 
accuracy  and  needs  a fairly  exacting  descriptive  technique. 

To  achieve  this,  direct  description  has  some  advantages  over 

indirect  description.  These  can  be  seen  in  a comparison  of  the 

two  approaches.  y,  S.  DEPT.  OF  hGRICOLTURc 

fiWlONAL  AGRICULTURAL  LIBRARY 

APR  3 01965 

C & R.KREP. 

1/  Statistician,  Central  States  Fordst  Experiment  Station, 

U.S.  Department  of  Agriculture,  Forest  Service. 


Indirect  Approach  that  of  selecting  a form  involving  a 


limited  array  of  commonly  used  transformations  of  the  independent 
variables  --- 

e.g.  Y = a + b^X  + b^X"^  + b^X^  + b^  logX 

fitting  this  form  to  a relevant  data  set  ---  statistically 

screening  components  for  those  most  able  to  reduce  the  variance 

about  regression  and,  hoping  the  remaining  ones  will  comprise 

2/ 

a form  reasonably  close  to  the  graphic  standard,— 

In  this  case,  useful  advance  knowledge  regarding  the  form  of 
the  relation  is  ignored.  The  amount  of  information  salvaged  from 
the  data  alone  depends  on  data  strength  and  on  the  degree  to  which 
the  limited  set  of  monomials  above  can  be  "molded"  into  a form 
like  the  standard. 

Direct  Approach  that  of  selecting  from  a large  array  of 

transformations  and  associated  coefficients  those  that  jointly  are 

3/ 

known  to  duplicate  the  graphic  form  with  reasonable  accuracy— 

fitting  the  resulting  algebraic  form  to  a relevant  data  set  by 


2_/  Where  the  graphic  form  is  developed  from  the  same  data  to 
which  it  is  being  fitted  by  Least  Squares,  the  propriety  of  statis- 
tical evaluation  may  be  questioned. 

The  standard  of  accuracy  is  set  by  the  analyst. 
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Least  Squares  and,  where  the  graphic  form  is  developed 

independently  of  the  data  to  which  it  is  fitted,  statistically 

2 

evaluating  the  results assuming  e^NID  (0,a  ). 

In  this  case  advance  knowledge  regarding  the  form  of  the 

relation  is  as  fully  utilized  as  is  practicable  in  establishing 

the  algebraic  form  so  that  fitting  the  algebraic  form  to 

relevant  data  will  tend  to  maximize  the  information  extracted. 

The  iterative  procedures  presented  in  this  paper  represent 

at  least  a rudimentary  effort  to  emulate  the  Direct  Approach  and 

4/ 

to  benefit  from  its  attributes.— 

It  should  be  noted  that  neither  this  nor  any  other  process 
offers  a panacea  for  curve  form  description  problems.  There  are 
likely  to  be  occasions  where  use  of  a Direct  approach,  as  demon- 
strated, will  yield  substantial  improvements  over  the  Indirect 

Approach  and  others  where  it  will  not.  Too,  success  may 

hinge  largely  on  the  ingenuity  of  the  user. 
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4/  Texts  often  contain  an  array  of  standard  algebraic  forms 
from  which  the  analyst  may  be  able  to  match  his  "expected" 
form  (1).  As  in  this  paper,  a more  basic  understanding  of  scaling 
and  of  algebraic  form  components  is  occasionally  sought  (2). 
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Examples  are  given  for  16  types  of  2 -dimensional  forms, 
two  3-dimensional  forms,  and  one  4-dimensional  form.  Brevity 
has  been  sacrificed  in  order  to  present  a uniform  approach  to  the 
description  problem.  But  given  understanding  of  the  basic  approach, 
you  will  be  in  a position  to  innovate  ---  the  field  is  virtually 
limitless . 

Relatively  high  powers  of  variables  are  used  on  occasion  to 
describe  sharp  curvature,  but  we  are  not  limited  here  except  by 
convention  and  by  computational  inefficiency.—^ 

Where  the  form  developed  involves  more  than  one  regression 
component,  be  sure  that  the  entire  function  will  be  suitably  smooth 
between  representative  control  points  from  the  model. 

Algebraic  forms  of  the  examples  have  been  left  unsimplified 
to  facilitate  understanding  of  the  source  of  scaling  coefficients. 
Simplification,  if  desired,  is  left  to  the  reader. 

The  more  complex  and  sometimes  more  flexible  descriptive 
system  involving  the  ratios  of  continuous  functions  should  be 
recognized  as  extremely  important,  although  not  covered  here. 
Examples  of  "ratio-forms"  are  given  by  Elderton  (3)  in  his  discus- 
sion of  Karl  Pearson's  basic  frequency  distributions. 


An  invaluable  aid  for  desk  computation  is  "Barlow's 
Tables  of  Squares,  Cubes,  Square  Roots,  Cube  Roots,  etc."  or 
some  similar  compilation. 
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THE  EXPECTED  FORM 


It  is  assumed  itl  this  paper  that  a graphic  form  has  already 
been  established  for  the  relation  involved.  EXTREME  CARE  is  urged 
in  developing  such  a form.  Before  looking  at  new  data,  construct 
an  expected  response  form  covering  extremes  of  the  independent 
variable(s).  Usually  you  can  be  quite  certain  as  to  the  general 
slope  and  value  of  the  curve  at  such  extremes,  e.g.,  flat,  tending 
to  zero  or  vertical,  tending  to  infinity.  Then,  what  happens  to 

the  curve  between  the  extremes  is  rather  strongly  suggested  

see  the  dotted  lines  in  the  graphed  examples  below  for  the  shape 
of  the  curve  between  the  extremes.  Points  in  the  range  of  the 


0 HI  0 Ml 

-X-  -x- 


independent  variable  at  which  important  changes  in  slope  occur  are 
sometimes  predictable  from  past  experience.  Exert  every  effort  to 
bring  known  information  to  bear  on  this  estimate  and  quantify  the 
relation  where  possible. 

Then  observe  the  shape  of  the  expected  curve  in  the  range  of 
the  independent  variate  that  is  relevant  to  the  new  data.  Describe 
this  curve  segment  algebraically  and  fit  it  through  the  new  data  by 
Least  Squares.  Evaluate  statistically  as  desired. 
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3-DlMENSIONS 


GRAPHIC  DEVELOPMENT  OF  SURFACES  p.  35 

ALGEBRAIC  DESCRIPTION: 


7 


TWO  DIMENSIONS 


TYPE  1 

INTRODUCTION  TO  SCALING 

Examples  of  positive  and  negative  slopes  from  various 
intercepts—^  are  presented  below.  Note  in  each  case  that  Y must 
be  described  in  terms  of  X,  Y-values  being  obtained  by  scaling  a 
given  value  of  X to  a corresponding  value  of  Y. 


Y=  A + B X 


Intercept  Y-value  at  X = 0 
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TYPE  2 


(oriented  at  Y = 0) 

We  might  start  here  by  seeing  how  well  some  scaled  powers  of 

X will  match  the  curve  or  model  at  representative  points  (dots  in 

2 

the  graphs  below).  For  example,  let's  try  X 

Before  scaling  After  scaling 


Scale  computation 

2 

Where  b = Y/X  , at 

7 / 

the  largest  value  of  X.— 


b = (12/16)  = .75 
and  Y = a + bX^ 

Y = 0 + .75  X^ 


l_l  Generally  a good  working  rule.  Although  sometimes  useful, 
scaling  at  some  other  value  of  X will  not  be  demonstrated  until 
later. 
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After  scaling  X , we  can  see  that  this  curve  is  too  flat  to 
represent  the  model  satisfactorily.  Because  higher  powers  of  X have 
more  curvature,  X and  X might  be  tried  next. 


Scale  computations 


After  scaling 


4 

For  X , 

b = (12/256)  = 

.0469 

For  X^, 

b = (12/4096)  = 

.00293 

And , 

Y = 0 + 

.0469 

Y = 0 + 

.00293  X^ 

(X) 


We  might  be  satisfied  with  the  proximity  of  Y = 0 + .0469  X 
to  the  model. 

But,  the  above  example  only  illustrates  the  relative  curvature 
of  some  integral  powers  of  X.  On  occasion,  it  may  be  found  that 
fractional  powers  of  X are  more  appropriate. 

8_/  Fractional  powers  of  X may  be  computed  from  X*^  = antilog  of  (n  log  X). 
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We  can  also  resort  to  additive  combinations  of  various  powers 


of  X.  Suppose  we  were  not  satisfied  with  Y = 0 + (.0469)  X^,  as  a 
description  of  the  original  model.  We  might  try  to  improve  that 
description  as  follows: 

From  an  examination  of  the  model,  it  appears  that  Y-values  for 
X = 0 through  X = 1.5  would  be  fairly  well  represented  by  a 
straight  line,  Y = a + bX  --- 
Where , 

b^  = (Y/X)  at  X = 1.5 
b^  = (.55/1.5)  = .367 
And , 

Y^  = 0 + .367  X 
Plotting  this  form  over  that 
of  the  model , (graph  at 
right)  we  see  that  the  quan- 
tities that  must  still  be 

added  to  match  the  model  are 
of  the  Type  I curve  form. 

So,  let's  add  a higher 
power  of  X as  before,  but 

scale  it  to  the  largest  additional  quantity  needed,  10.53  here. 
Also,  we  must  refrain  from  adding  large  quantities  to  the  already 

satisfactory  straight  line  to  the  left  of  X = 1.5  . Then 

the  power  of  X we  decide  to  use  must  be  associated  with  a small 
"b"  to  suppress  scaled  values  at  X = 1.5  and  below. 
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Let's  try  and  — For  X^,  = (10.53/256)  = .04113 

— For  X^,  = (10.53/4096)  = .002571 

and  = 0 + .04113  X^ 

= 0 + .002571  X^ 

Adding  these  forms  to  the  straight  line  already  plotted,  we 
have 


(X) 
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From  the  graph  above,  we  can  readily  see  that  in  matching  the 

model,  the  better  of  the  two  curves  is  = .367  X + .002571 

2 4 6 

and  that  it  is,  in  fact,  better  than  Y = X^,  X , or  X alone 
(tried  previously). 

A bookkeeping  system  follows  for  the  computations  - — 


Model 

Values 

Added  amt. 
to  match 
model 

Difference 
from  model 

X 

Y 

(.55/1.5)  X 

= .367  X 

Y-Yi 

x" 

(10.53/4096)  X^ 
= .002571  X^ 

+ 

>* 

1 

0.0 

0.00 

0.00 

0.00 

0 

0.00* 

0.00 

1.0 

0.30 

0.37 

-0.07 

1 

0.00* 

-0.07 

1.5 

0.55 

0.55 

0.00 

11 

0.03* 

-0.03 

2.0 

1.00 

0.73 

0.27 

64 

0.16 

0.11 

3.0 

3.00 

1.10 

1.90 

729 

1.87 

0.03 

4.0 

12.00 

1.47 

10.53 

4096 

10.53 

0.00 

* negligible  additions  to  for  X = 1.5  or  less 


Note  that  = 0.367  X + 0.002571  X^  differs  from  the  model 


by  0.11  or  less  for  the  X-values  examined  over  the  range  of  X. 
Our  final  curve  form  here  is  Y = a + b^X  + b^X^.  A more  exact 
description  of  the  model  might  or  might  not  be  found  by  trying 
other  forms.— ^ 


£/  In  place  of  0.367  X,  we 
power  curve , such  as , b^X^ • ^ 


might  start  with  a fairly  flat 
perhaps  next  adding  b^X^  etc. 
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TYPE  2 


(not  oriented  at  X,Y  = 0) 

Where  the  model  intersects  the  Y-axis  at  some  point  other 
than  zero,  the  curve  development  procedure  would  be  identical  to 
that  of  the  foregoing  case,  but  straight-line  and  curved  forms 
would  now  be  added  to  the  intercept  value,  Y = a. 


(X) 


As  before,  final  form  adopted  here  is  Y = 


a + b^X  + b^X 


14 


TYPE  3 


The  procedure  here  is  basically  the  same  as  that  outlined 
for  T3rpe  2,  Note  the  transformation  of  X,  (4-X) , for  the  last 
two  terms  — =■  to  make  additions  to  regression  larger  as  X approaches 
zero.  About  the  same  effect  could  be  obtained  using  l/x’^,  X > 0, 
n > 1,  but  this  transformation  is  a little  more  difficult  to 
compute.  Also,  starting  the  building  process  at  a = 14  would 

permit  use  of  the  untransformed  ‘'X"  in  Y = 14+(-b^)X+b2X” 

but  text  continuity  is  best  served  here  and  subsequently  by 
adopting  (4-X)  when  appropriate  either  transform  is  fully 
acceptable. 


(X) 

4 3 2 10 

(4-X) 


Y=2  + .367{4-X) 
+.00  2 57/(4-X) 
Y = 2 + .36  7 (4-X) 
Y=2 


6 


The  final  form  description  adopted  in  the  above  example  is 
Y = a + b^(4-X)  + b2(4-X)^  0 ^ X < 4 
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TYPE  4 


The  procedure  here  is  basically  the  same  as  that  outlined  for 
TYPE  2 except  that  we  are  working  from  high  to  low  values  of  Y. 

Note  the  transformation  of  X,  (4-X) , for  the  last  two  terms  to 
make  reductions  from  the  constant  (12)  larger  as  X approaches  zero. 

Y = I2 


(X) 

4 3 2 10 

(4-X) 


The  final  curve  fom  description  adopted  in  the  above  example 
is  Y = a + b^(4-X)  + b2(4-X)^,  0 < X ^ 4 

Another  possibility  here  is  to  start  with  Y = 0,  adding  scaled, 

positive  fractional  powers  of  X to  arrive  at  a form 

like  that  of  the  model.  But  once  again,  X^^’^  is  more  difficult  to 
compute. 
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TYPE  5 


The  procedure  here  is  basically  the  same  as  that  outlined  for 
TYPE  2,  except  that  we  are  working  from  high  to  low  values  of  Y. 


(X) 


The  final  curve  form  description  adopted  in  the  above  example 
is  Y = a + b^X  + b^X^ 
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TYPE  6 


Consider  this  form 
in  two  pieces. 


1.  Describe  part  (A)  as  per  TYPE  4. 

2.  Subtract  a high  power  of  X from  the  (A)  form  to  give  the 

sharp  downward  trend  of  part  (B)  

For  example,  using  the  model  below: 


A B 


* 


4 3 2 10 

(4-X) 


* From  extension  of  the  model  trend. 
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To  match  part  (A)  of  the  original  model,  we  transform  X to 
2 

(4-X)  and  scale  this  quantity  to  10.5  (12.5  - 2.0)  at  the 

2 

largest  value  of  (4-X).  We  find  that  Y = a - b(4-X)  is  fairly 
suitable.  Scaled,  it  is  Y = 12.5  + (-10.5/16)  (4-X)^. 

To  match  part  (B)  of  the  model,  we  must  find  some  high  power 
of  X that,  when  scaled  to  8.5  and  subtracted  from  our  description 

for  part  (A),  will  subtract  little  or  nothing  from  X=0  through  X=3. 

20  9 -11 

Try  X , then  b = (-8.5/1100  x 10  ) = -0.773  (10  ) and  the 

complete  model  can  be  described  as  --- 

2 -1120 

Y = 12.5  - 0.656(4-X)  - 0.773(10  )X  -—  0 < X < 4 


19 


Summary  of  Computations 


Part  (A) 


Model 

^"2 

X Y 

a* 

Y 

- ’'i 

(4-X) 

(4-X)^ 

-10.5 

16 

(4-X)^ 

Y-CY^+Y^ 

0 2.0 

12.5 

• 

10.5 

4 

16 

10.50 

0 

1 6.8 

12.5 

- 

5.7 

3 

9 

- 

5.90 

.20 

2 10.2 

12.5 

- 

2.3 

2 

4 

- 

2.62 

.32 

3 11.8 

12.5 

- 

.7 

1 

1 

- 

. 66 

- .04 

4 4.0 

12.5 

- 

8.5 

0 

0 

- 

0 

- 8.50 

Part  (B) 


Model 

diff.  from  model 

Y V 

-8.5  20 

Y-(Y^+Y2+Y3) 

A 1 

I.IOOMMM 

0 2.0 

0 

0 

0 

1 6o8 

1 

o 

0 

.20 

2 10.2 

1,049 

(lop 

0 

.32 

3 11.8 

3,487 

(lob 

- .03 

- .01 

4 4.0 

1,100 

aob 

- 8.50 

0 

* From  extension  of  model  trend 


Providing  the  above  "difference  from  model" 
excessive,  the  accepted  final  form  would  then  be 


is  not  considered 


Y = a + b^(4-X)^  + b^(X)^^ 


---  0 < X < 4 
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TYPE  7 


The  procedure  outlined  for  TYPE  6 is  almost  identical  to  the 
one  applied  here  — - follow  the  steps  in  the  graph  below. 


(X) 

y.  * ® *FROM  EXTENSION  OF 

^ MODEL  TREND 


Then,  the  final  form  here  is  --- 

Y = a + + b2(4-X)^®  0 <X  < 4 
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TYPE  8 


The  procedure  outlined  for  TYPE  6 is  very  similar  to  the  one 
applied  here  follow  the  steps  in  the  graph  below. 


14 
12 
10 
8 

(Y)  6 

4 
2 


(X) 

4 3 2 10 

{4-X) 

* FROM  EXT  ENSIGN  OF 
MODEL  TREND 

The  final  form  is  then  

2 20 

Y = a + b^(4-X)  + b^X  — 0 < X <4 
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TYPE  9 


The  procedure  outlined  for  TYPE  6 is  very  similar  to  the  one 
applied  here  ---  follow  the  steps  in  the  graph  below. 


(X) 

4 3 2 10 

{4-X) 

* FROM  EXT  EN  SION  OF  MODE  L TREND 


The  final  form  adopted  is  then  

2 20 

Y = a + b^X  + b2(4-X)  ---  0 < X < 4 
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TYPE  10 


Where  the  curvature  on  either  side  of  the  minimum  is  more  or 
less  symmetrical,  apply  TYPE  2 processes  to  the  transformation  of 
X given  below. 

Determine  *'X*'  at  the  minimum  value  of  "Y"  and  subtract  from 

10/ 

all  X-values.  Then,  using  the  absolute  values — of  these  differ 
ences,  proceed  to  locate  descriptive  terms  for  either  side  of  the 
model . 

For  example: 


-3-2-10  I 2 3 
(X-3) 


10/  Use  of  the  absolute  value  permits  adoption  of  any 
integral  or  fractional  power  of  the  transform  through  elimination 
of  negative  algebraic  signs. 
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Summary  of  Computations; 


Model 

^^2 

X 

Y 

a* 

(Y-Y^ 

(X-3) 

|x-3|2 

0 

12.0 

2 

10.0 

-3 

9 

6.98 

1 

5.1 

2 

3.1 

-2 

4 

3.10 

2 

2.8 

2 

.8 

-1 

1 

.78 

3 

2.0 

2 

.0 

0 

0 

. 00 

4 

2.8 

2 

.8 

1 

1 

.78 

5 

5.1 

2 

3.1 

2 

4 

3.10 

* The  minimal  point  of  the  graphed  curve 


Model 

’^3 

diff.  from  model 

X 

Y-CY^+Y^) 

X 

1 

u> 

00 

3.02.  .18 
6561 

Y-CY^+Y^+Y^) 

0 

3.02 

6561 

3.02 

0.00 

1 

.00 

256 

.12 

-.12 

2 

.02 

1 

.00 

.02 

3 

.00 

0 

.00 

. 00 

4 

.02 

1 

.00 

.02 

5 

.00 

256 

.12 

-.12 

The  final  form  adopted  is 

Y = a + bjx-sl^  + b^jx-sl  ^ ---  0 < X 6 


25 


TYPE  11 


Where  the  curvature  on  either  side  of  the  maximum  is  more  or 
less  symmetrical,  apply  TYPE  2 processes  to  the  transformation  of 
X given  below. 

Determine  ’’X"  at  the  maximum  value  of  ”Y*'  and  subtract  from 
all  X-values,  Then,  using  the  absolute  values  of  the  differences, 
proceed  to  locate  descriptive  terms  for  either  side  of  the  model 
as  shown  for  TYPE  2.  Note  that  we  are  working  down  from  the  inter- 
cept, Y = 10,  here.  The  computations  are  directly  analagous  to 
those  for  TYPE  10, 


12 


0 12  3 4 5 6 
(X) 

-3-2-1  0 I 2 3 
(X-3> 


6 


The  final  form  adopted  is 


Y = a + b^  iX-3l^  + b^  lX-3l®  — 0 IX  < 6 
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TYPE  12 


When  trends  on  opposite  sides  of  the  inflection  point,  I— ^ , 

are  reasonably  symmetric  in  reverse  directions,  the  following 
procedure  may  lead  to  a fairly  good  description  of  the  model. 

1.  Let  a straight  line  represent  as  much  of  the  central 
section  of  the  model  as  seems  appropriate.  Find  graphically  the 
Y-value  of  this  line  at  X = 0 and  also  the  slope  of  the  line. 

See  graph  (A)  in  the  example  that  follows. 

2.  It  will  be  seen  from  graph  (A)  that  it  is  necessary  to 
add  a curvilinear  quantity  to  the  lower  end  of  the  line  and  to 
subtract  a similar  amount  from  the  upper  end  of  the  line  to  match 
the  model . 

Transforming  "X"  to  X - (X  at  the  inflection  point),  we  have 
a series  of  numbers  that  are  equal  but  opposite  in  sign  on  each 

side  of  the  inflection  point  X^.  The  sign  differences  are  main- 
tained whenever  the  integral  power  of  (X=X^)  is  odd  

e.g.:  (X-X^)\  (X-X^)^,  (X-X^)^ 

or  when  the  denominators  of  fractional  powers  are  odd  integrals  (4) 
e.g.:  (X-X^)^/^,  (X-X^)^^/^, 


11/  Note  inflection  point,  I,  in  graph  (A)  which  follows. 
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to  add  and  subtract 


Then  we  can  use  the  above  powers  of  (X-X^) 

from  opposite  ends  of  the  straight  line,  as  indicated  in  graph  (B) 

n 

below.  Note  that  the  scaling  coefficient  for  (X-Xj.)  must  be 
small  enough  to  suppress  additions  to  the  straight  line  values 
that  already  represent  the  model  with  suitable  accuracy. 


(X-3) 
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Summary  of  Computations ; 


Model 

Y 

1 

^^2 

^3 

diff . fr. model 

X 

Y 

13.2*„ 

Y-CY^+Y^) 

(X-3) 

(X-3)^ 

(2.6)(X-3)'^ 

Y-(Y^+Y2+Y^) 

6 * 

-27 

0 

2.0 

-0.6 

0.  0 

2.60 

- 3 

-27 

2.60 

0.00 

1 

2.4 

- .6 

2.2 

.80 

- 2 

- 8 

.77 

.03 

2 

3.8 

- .6 

4.4 

.00 

- 1 

- 1 

.10 

- .10 

3 

6.0 

- . 6 

6.6 

.00 

0 

0 

.00 

.00 

4 

8.2 

- .6 

8.8 

.00 

1 

1 

- .10 

.10 

5 

9.6 

- .6 

11.0 

- .80 

2 

8 

- .77 

- .03 

6 

10.0 

- .6 

13.2 

-2.60 

3 

27 

-2.60 

.00 

* determined  graphically 


The  final  form  adopted  is 

Y = a + b^X  + b2(X-3)^  0 < X < 6 
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TYPE  13 


The  procedures  for  TYPE  12  are  directly  applicable  here, 
excepting  that  the  sign  of  the  straight  line  becomes  negative. 
For  example: 


Y=  Q + b|X  -»*b2(X~3)^ 


(X) 

-3-2-10  1 2 3 

(X-3) 


Summary  of  Computations: 


Model 

^"2 

^3 

o 

X Y 

a* 

[^] 

* 

X 

Y-(V^2) 

(X-3) 

(X-3)^ 

2 (X-3)'" 
“27 

Y-(Y^+Y2+Y^) 

0 14.0 

12 

0. 

00 

2.00 

-3 

-27 

2.00 

0.00 

1 11.1 

12 

-1. 

33 

.43 

-2 

- 8 

.59 

-.16 

2 9.3 

12 

-2. 

67 

-.03 

-1 

- 1 

.07 

-.10 

3 8.0 

12 

-4. 

00 

. 00 

0 

0 

. 00 

.00 

4 6.7 

12 

-5. 

33 

.03 

1 

1 

- .07 

.10 

5 4.9 

12 

-6. 

67 

-.43 

2 

8 

- .59 

.16 

6 2.0 

12 

-8. 

00 

-2.00 

3 

27 

-2.00 

,00 

* determined  graphically 

The  final  form  adopted  is 

Y = a + b^X  + b^CX-S)^  ---  0<  X <6 
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TYPE  14 


The  procedures  for  TYPE  13  are  directly  applicable  here 
excepting  that  the  sign  of  in  b^CX-X^)^  is  positive. 

TYPE  15 


The  procedures  for  TYPE  12  are  directly  applicable  here 
excepting  that  the  sign  of  b2  in  b2(X-X^)*^  is  positive. 
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DRAINAGE  RATE.cc/SEC. 


TYPE  16 


The  model  below  is  a'  mixture  of  TYPES  3 and  7 (pages  14  and 
20  respectively)  applied  to  a particular  research  generated  graphic 
form. 


A tabular  summary  of  the  algebraic  form  development  follows. 
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SUMMARY  OF  COMPUTATIONS 
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NOTES 


1.  Here,  reciprocals  of  the  independent  variable,  t,  replace  powers 

of  the  linear  transposition  (4-X)  shown  for  TYPES  3 and  7. 

The  reciprocals  happened  to  be  more  accurate  in  this  case. 

2.  describes  the  back-slope  of  the  model  and  is  scaled  to  D +.35 

-4 

at  t = 10.00  ---  scaling  coefficient  = (9.70  +.35)/1.000  X 10  = 

4 

10.05  xlO  . The  .35  was  added  at  this  point  because  a trial  scaling 
using  D alone  showed  the  right  tail  to  be  about  .35  cc/sec.  high. 


3.  In  X^  (or,  X^-.35)  the  -.35  corrects  for  the  high  right 

tail  of  the  back-slope. 

2 

4.  And  last,  appropriately  increasing  amounts,  -.549/(t-9)  , were 

subtracted  from  X^  as  t became  smaller.  "9"  was  used  as  the  point  of 

2 

origin  for  the  t-transform,  (t-9)  , because  9 was  the  origin  of  the 
model.  The  scaling  coefficient  for  X^  in  this  case  was  based  on  the 

difference  between  X2  and  D of  the  model  at  t = 9.67  X^  = 11.15, 

D = 9.90,  X^-D  = 1.25  ---  l/(t-9)^  = 2.25  ---  coeff.  = 1.25/2.25  = .551 

finally  adjusted  to  .549  for  slight  benefit  in  the  region  to  the 

left  of  maximum  D. 
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THREE  DIMENSIONS 


GRAPHIC  DEVELOPMENT  OF  SURFACES 

Before  discussing  the  mathematical  description  of  multi- 
dimensional models,  it  will  probably  be  desirable  to  consider  the 
graphic  development  of  surfaces  from  sets  of  data.  A fairly 
straight-forward  procedure  is  presented  in  the  Appendix,  beginning 
on  page  1. 

ALGEBRAIC  DESCRIPTION 

The  analyst's  objectives  and  the  degree  of  his  confidence  in 
the  adopted  model  will  probably  control  the  extent  of  his  algebraic 
description  effort.  Just  as  for  2-dimenSional  models,  we  will 
consider  the  case  where  a fairly  accurate  description  of  a well- 
established  model  is  desired  and  where  complexity  of  algebraic 
description  is  not  a deterrent  in  either  fitting  or  use  of  the  model. 

Example  #1.  It  will  be  seen  here  that  the  algebraic  form 
derived  is  uncommonly  complex,  involving  five  transformations  of 
the  two  independent  variables,  which  is  the  price,  in  this  case, 
for  a reasonable  facsimile  of  the  model.  Given  additional  effort, 
it  is  possible  that  a simpler,  more  accurate  form  could  be  developed. 
But  the  example  is  retained  because  it  is  fully  applicable  "as  is" 
and  demonstrates  the  descriptive  process  well. 

The  solid-line  figure  in  the  following  graph  is  the  model  to 
be  described. 
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MODEL  SURFACE 
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MODEL  SURFACE  VS.  A DERIVED  MATHEMATICAL 
DESCRIPTION  OF  THAT  SURFACE 


To  describe  the  model  algebraically,  start  by  selecting  enough 
surface  values  over  the  ranges  of  both  independent  variates  (X,Y) 
to  specify  major  surface  trends  adequately.  But  keep  the  number  of 

such  values  to  a minimum  for  twelve  (X,Y) -combinations  were 

adopted  here. 

Then  proceed  to  find  form  components  that  jointly  describe  an 
acceptably  smooth  surface  passing  through  the  Z-points  adopted.  The 
final  components  are  shown  below  in  tabular  form  in  preference  to 
lengthy  verbal  description.  Note  how  successive  components  describe 
additional  portions  of  the  model  much  as  they  did  for  2-dimensional 
cases.  Components  parallel  those  for  curve  TYPES  3 and  9.  Other 
surfaces  may  involve  any  or  all  of  TYPES  1-16. 

Two  added  features:  scaled  transformations  of  the  Y-values  are 
used  to  suppress  additions  to  portions  of  the  surface  for  which  satis- 

12/ 

factory  X-descriptive  terms  have  already  been  located — and 

« 

(X,Y)  -cross  products  are  developed  to  describe  joint  changes  in 
X and  Y. 

Each  operation  must  be  carefully  thought  through  in  order  to 
follow  the  steps  outlined.  The  exponent  adopted  for  any  term  is 
the  best  of  perhaps  five  or  six  tried.  The  entire  description 
here  required  two  days  for  development. 

A graphic  portrayal  of  the  regression  components  follows  the 
derivation  tables. 


12/  Both  X and  Y are  independent  variates  here. 
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DEVELOPMENT  OF  ALGEBRAIC  DESCRIPTION 
OF  SURFACE  MODEL 


Model 

A 

„1.85 

r .89 

„1 . 85 

Y 

Y X 

X 

[36.6 

X 

15 

0 1 

0.000 

1.00 

0. 

024 

0 

3 

.003 

7.63 

* 

186 

0 

5 

.019 

19.64 

* 

478 

0 

7 

.000 

36.60 

• 

890 

0 

15  1 

.410 

1.00 

024 

1 

3 

.183 

7.63 

* 

186 

1 

5 

.479 

3 

' 19.64 

478 

1 

7 

.890 

36.60 

, 

890 

1 

25  1 

1.200 

1.00 

024 

1.667 

3 

.403 

7.63 

, 

186 

1.667 

5 

1.049 

19.64 

478 

1.667 

7 

2.340 

36.60 

• 

890 

1.667 

(7-X) 


10 


.386 


0,466M 


(7-X) 


10 


15 


0 

0 

0 

0 

0.024 

.186 

.478 

.890 

.040 

.310 

.797 

1.484 


0.000 

.003 

.019 

0 

.386 
■ .003 
.001 
0 

1.160 

.093 

.252 

.856 


60 ,466M 
1,049M 
1,024 
0 

60,466M 

1 , 04 9M 
1,024 
0 

60,466M 

1,049M 

1,024 

0 


0 

0 

0 

0 

.386 

.007 

0 

0 

.643 

.011 

0 

0 


A constant  the  lowest  value  for  the  surface. 

7^1  Curve  Type  3 ---  well  described  by 

3^/  Curve  Type  9 described  by  ---  note  that  Y is  used  to  suppress 

additions  to  + Z^,  already  satisfactory  at  Y = 0. 
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,3.6 


.856 


1102.5 


,3.6 


■ y“ 

7 

r ^ 1 

"y 

15 

L.  — 

35.78) 

15 

BC 


15 


25 


0.000 

.003 

.019 

.000 

.000 

-.010 

.001 

.000 

.517 
.082 
.252  1/ 
.856 


1.0 

52.2 

328.3 

1102.5 

1.0 

52.2 

328.3 

1102.5 

1.0 

52.2 

328.3 

1102.5 


0.001 

.041 

.255 

.856 

.001 

.041 

.255 

.856 

.001 

.041 

.255 

.856 


0 

0 

0 

0 

1 

1 

1 

1 

35.7  8 
35.78 
35.7  8 
35.78 


0 

0 

0 

0 

.028 

.028 

.028 

.028 

1 

1 

1 

1 


0 

0 

0 

0 

0 

,001 

.007 

.024 

.001 

.041 

.255 

.856 


Y X 

^6 

2/ 

(7-X)^° 

.516“ 

(7-X)^°C 

60,466M 

0 1 

0.000 

60,466M 

0 

0.000 

3 

.003 

1,049M 

0 

.003 

5 

.019 

1,024 

0 

.019 

7 

.000 

0 

0 

.000 

15  1 

.000 

60,466M 

.014 

-.014 

3 

-.011 

1,049M 

0 

-.011 

5 

-.006 

1,024 

0 

-.006 

7 

-.024 

0 

0 

-.024 

25  1 

.516 

60 ,466M 

.516 

.000 

3 

.041 

1,049M 

.009 

.032 

5 

-.003 

1,024 

0 

-.003 

7 

.000 

0 

0 

.000 

1^/  Curve  Type  9 described  by  -- 

suppress  additions  to  Z.  + Z„  + Z.  + Z,, 
Y = 0 and  Y = 15. 


note  that 


-_f]7 
15 


is  used  to 


already  satisfactory  at 


7^1  These  are  the  aggregate  differences  between  the  model  values  and  the 

mathematical  description  note  the  small  differences  in 

comparison  to  the  original,  corresponding  model  Z values. 
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FINAL  ALGEBRAIC  DESCRIPTION 


Z = .51  + .00896(7-X)^*^  + 


.00077 


.3.6 


j-— 


-j.  |.0007764(X) 

.027988115J  1— O^Xi7 


. 02432 (X)^*®^  [l5 
.027988  {is]  ] + 


.6384(10'®)(7-X)^® 


.8534(10'®)(7-X)^® 


This  form  is  graphed  over  that  of  the  model  in  dotted  lines 

for  visual  comparison see  the  foregoing  graph. 

Also,  see  the  following  graph  for  contribution  to  regression 
by  components  of  the  form  above. 
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HAUL  TIME  CMIN.) 

PER. MBM  PER.  too  FOOT 


CONTRIBUTION  TO  REGRESSION 
BY  COMPONENTS  OF 
THE  DERIVED  FORM 


LOAD,  100  BO. FT. 
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Example  #2.  The  model  that  follows  was  developed  from  a very 
large  number  of  observations.  Feeling  that  this  model  is  likely 
to  be  very  representative  of  the  population  form,  one  might  adopt 
it  for  testing  in  a new  sample  from  a similar  population. 

Again,  with  the  major  objective  of  developing  a reasonably 

13/ 

accurate  algebraic  description,  our  computations — follow  about 
the  same  pattern  as  previously,  so  some  of  the  detail  of  scaling 
is  omitted.  Instead,  a brief  verbal  description  of  the  individual 
algebraic  steps  and  graphs  showing  regression  components  are  included. 

Note  that  "absolute  value"  is  used  to  obtain  increases  in  the 
same  direction  on  either  side  of  the  curve  of  minimums  (trough  base) 

see  "C" . Also,  scaling  of  a transformation  may  be  based  on 

slightly  less  or  more  than  maximum  values,  so  that  an  average  coeffi- 
cient may  be  used  to  represent  two  or  more  groups  of  representative 
points  of  the  model.  This  is  true  in  the  scaling  of  both  "C"  and 
"G"  in  the  computations  that  follow. 


13/  Computations  follow  the  model. 
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I.  - REPLOTTEO 


HEIGHT 


H.  BASIC  TROUGH, SCALED 


HEIGHT 
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SPECIFIC  GRAVITY 


SIDE  VIEW 


TIT,  ELEVATION  OF  TROUGH 

ABOVE  CONSTANT,  .311 


.08i 


.06- 

.05- 

.04- 

.03- 

.02- 

.01- 

0- 


1.6  X I0”^{|a-[25  + (-^— )]jH  {50  — H) 
■+•.43x10“'®  (50-H)® 


HEIGHT 
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ADDING  THE  C ONSTANT,  . 3 1 1 , TO  THE  ABOVE  FORM  BRINGS 
THE  DERIVED  FORM  REASONABLY  CLOSE  TO  THAT  OF 
THE  ORIGINAL  MODEI SIMPLY  SHOWN  IN  THE 

SIDE  VIEW  below: 
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FOUR  DIMENSIONS 


A fairly  simple  case  will  be  presented  here  to  demonstrate 

application  of  the  description  system  to  representative  points 

14/ 

from  three  model  surfaces. — Note  that  each  surface  is  a point 
in  the  fourth  dimension.  The  continuity  between  surfaces  will  be 
provided  by  the  term  or  terms  in  the  algebraic  description  that 
involve  the  fourth  variable.  As  before,  take  care  that  the  terms 
selected  do  not  move  unrealistically  between  the  representative 
points  to  which  you  scale. 

The  case  at  hand  is  based  on  a large  amount  of  actual  data 
from  forest  fire  research.  Surfaces  were  developed  from  data  by 
procedures  shown  in  the  appendix.  An  additional  refinement,  that 
of  smoothing  number  of  fires  for  identical  (SMC-BUI)  points  from 
each  surface,  was  employed  to  arrive  at  adjusted  responses  for 
evenly  distributed  points  in  wind  velocity.  The  three  model  sur- 
faces below  were  the  result. 


14/  The  same  principles  can  be  applied  to  five  or  more 
dimensional  cases. 
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^ FIRES  PER  DAY 


3 


MODEL  DEVELOPED  GRAPHICALLY  FROM 
DATA  MEANS  FOR  GREEN 


SMC  JJ 


SMC 


SMC 


WIND  VELOCITY,  (mph) 


JJSTtCK  MOISTURE  C ON  T E N T , 
BUILD-UP  INDE  X 


An  algebraic  description  of  this  relation  is: 

# Fires  = . 282x10“^ + o477xl0“^(A^°  5 < C < 20 

shown  in  the  computations  which  follow., 

The  differences  between  representative  points  from  this  form 
and  corresponding  scaled  points  of  the  model,  are  not  discernible 
at  the  scale  shown. 


100  ^ 
BUI 
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DEVELOPMENT  OF  FORM  FOR  NUMBER  FIRES  PER  DAY 
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NOTES  ON  THE  FOREGOING  FORM  DEVELOPMENT: 


lo  selected  because  (D^*^  at  C = 5)  vs,  (D^*^  at  C = 10) 

has  about  the  same  ratio  or  slope,  as  at  (A  = 0,  B = 100, 
C = 5)  vs.  (A  = 0,  B = 100,  C = 10.)  e.g. 


i c 

3 1 

D ‘ ratio 

Yq  ratio 

5 

10 

4425\ 

\ 3.5 
1259X 

1.25v 

\3.6 

.35/ 

When  4425  is  scaled  to  1.25  — (1.25/4425)  4425  = (.282x10"^) 4425 
3 1 

= 1.25,  D * for  C = 10,  1259,  must  also  scale  close  to 

the  value  desired,  .35---(. 282x10  ^)  1259  = ,36 
3 1 

2.  D * is  multiplied  by  B in  forming  to: 

3.1 

(a)  suppress  D at  B = 0 

(b)  maintain  linear  additions  with  a change  in  B. 

3.  Y^  shows  the  amount  of  Y^  yet  to  be  added  in  the  algebraic 
description  of  the  model. 

4.  The  ratio  between  Y^  for  C = 5 and  C = 10  is  .90/. 17  = 5.3 

for  A = 6 and  2. 95/. 47  = 6.3  for  A = 12.  It  was  decided  here 

n 4.3 

to  adopt  an  average  ratio  for  D of  5.7,  provided  by  D at 

C = 5 and  10— (114,070/19,950)  = 5.7 
4 3 

Then,  D * scaled  to  .90,  will  yield  close  to  the  desired 
value,  .17,  for  C = 10  at  A = 6 — -and  scaled  to  2.95,  will 
yield  values  close  to  the  desired  value,  .47jfor  C = 10  at 
A = 12. 

5.  The  ratio  between  Y^  for  C = 5,  B = 100,  A = 6 and  12  is 

(2.95/. 90  = 3.3)  and  at  C = 10  is  (.47/. 17  = 2.8).  An  aver- 


50 


age  ratio  of  3.0  was  adopted  and  is  provided  by  A * 

(53.3/17.6  = 3.0). 

1 6 

Scaling  A ' at  (A  = 12,  C = 5)  will  then  yield  close  to  the 
desired  value,  .90  and  similarly  to  .47  and  .17,  at 


(A  = 6,  C = 5). 

X 6 3 

The  product  (A  * BD  ’ ) then  yields  a transformation  that 
suppresses  scaled  values  to  zero  at  the  appropriate  points  and 
allows  reasonable  additional  curvature  for  both  A and  D.  This 
composite  quantity  was  then  scaled  to  2.90  instead  of  2.95  in 
order  to  distribute  some  of  the  differences  from  more  evenly 
— - (2. 90/607, 993, 100)E  = ( .477xlO“®)E. 


7.  The  happenstance  presence  of  rather  uniform  slope  changes  in 
A = 6 and  12  permitted  the  ratio  averaging  technique  and 
resulted  in  a brief,  but  fairly  accurate  description  of  the 
model.  A more  exacting  description  could  probably  be  compiled 
with  the  use  of  separate,  algebraic  additions  for  each  change 
in  trend  encountered. 
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APPENDIX 


GRAPHIC  DEVELOPMENT  OF  SURFACES 
Any  relation  involving  two  independent  variables  can  be  shown 
as  a surface---more  than  two,  as  a series  of  surfaces---each  at 
specified  values  of  the  remaining  variables. 


For  example,  a relation  involving  the  dependent  variable  Z 
and  the  independent  variables  W,  X,  and  Y might  be  presented  as 
below. 


XXX 


Then  it  is  only  necessary  to  understand  the  derivation  of  a 
graph  for  a single  surface  in  order  to  make  the  above  presentation. 
Such  a derivation  scheme  follows  for  cases  where  there  is  a reason- 
ably wide  dispersion  of  observations  over  the  range  of  the  inde- 
pendent variables.  Without  such  dispersion,  other  techniques  are 
more  appropriate. 

Given:  Observations  of  Z on  W,  X and  Y 
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Procedure : 


(a)  Rate  the  independent  variables  as  to  suspected  degree  of 
importance  in  their  effects  on  Z. 

Separate  the  range  of  the  least  important  variable,  say 

''W" , into  four  parts thru  such  that  all  observations 

are  about  equally  distributed  among  them. 

(b)  Treat  each  W-group  as  follows: 

Separate  the  range  of  the  next  most  important  variable, 
say  "X",  into  four  parts  thru  X^  such  that  the  obser- 
vations are  about  equally  distributed  among  them  find  the  average 

X-value  for  each. 

(c)  For  each  X-group,  fit  a line  through  the  Z-values  plotted 
over  corresponding  Y-values , being  careful  to  adhere  to  expected 
trends . 


Repeat  for  remaining  X-groups. 


(d)  Transfer  these  curves  to  a 3-dimensional  graph.  Be  sure 
to  plot  each  curve  over  the  proper  average  X-value  obtained  in  (b). 


O 

Plot  the  3rd  dimension  at  a 45  angle  to  the  1st  dimension. 


(e)  Cross-curve  lightly,  as  below,  balancing  + and  - deviations 
as  nearly  as  is  practical. 


^ “ SMOOTHED  Z -VALUES 
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(f)  Using  the  cuin^ed  Z-values  above,  fit  a final  smooth 
curve  through  the  four  points  of  each  X-value  to  complete  the 
surface,  as  below; 


(g)  However,  it  may  be  desirable  to  obtain  surface  Z-values 
for  XY-values  other  than  those  represented  (or,  to  interpolate). 

It  is  then  necessary  to  locate  the  desired  intercepts  on  existing 
curves  and  to  join  these  intercepts  with  harmonizing  curves. 

The  surface  value  for  the  new  combination  of  X-Y  might  then  appear 
as  "a"  in  the  graph  below. 
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(h)  A further  refinement  can  be  applied  by  repeating  (a) 
through  (g)  switching  X and  Y wherever  either  or  both  appear  in 
the  directions.  You  then  have  derived  two  surfaces.  It  but 
remains  to  obtain  values  from  each  surface  for  each  of  say, 

sixteen  identical,  XY-values  and  to  average  each  pair  of 

corresponding  readings.  The  average  surface  may  then  be  plotted 
from  these  new  points. 

The  foregoing,  with  or  without  averaging  (as  time 
permits),  is  then  repeated  for  ^2  ^ and  4 

Final  surface  values  may  then  be  obtained  for  all  64 
WXY  combinations  and  a formula  derived  for  the  graphic  relation 
established. 
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